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Possible existence of the primordial magnetic fields has affected the structure formation of the 
universe. In this paper it is shown that the initial conditions for density perturbations with mag- 
netic fields derived in previous works are inconsistent with Einstein equations. We find that this 
inconsistency arises due to the unwanted cancellation of contributions from the magnetic fields and 
primordial radiations. A complete set of equations and consistent initial conditions in the long 
wavelength limit are given with an explicit derivation in the covariant approach with CDM frame, 
by newly taking into account a non-relativistic matter contribution in the radiation dominated 
era. By solving these equations numerically, we derive the angular spectrum of cosmic microwave 
background anisotropics and the matter power spectrum with magnetic fields. We find that the am- 
plitude of the angular power spectrum of CMB anisotropics can alter at most a order of magnitude 
at I < 4000 compared with the previous results in the literature. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 



There are mounting evidences that the large scale magnetic fields are present in various objects in the universe. 
Not only galaxies, but also clusters of galaxies contain their own magnetic fields with the field strength of ~ 10~^ 
Gauss and the coherence length of 1 — 10 kpc (for a review, Furthermore, there have been some observations 
which indicate that they exist even in larger scales, such as in superclusters 

Yet the origin of such large scale magnetic fields is still a matter of debate. Magnetic fields in spiral galaxies are 
assumed to be continuously generated and maintained by dynamo mechanism Q. However, one still needs to explain 
the origin of seed fields necessary for dynamo action to take place. Astrophysical origins of such seed fields, often 
involving stellar activities or the Biermann battery in non-adiabatic processes 0, S @ j may explain the strength and 
the total amount of magnetic fields with help from the dynamo mechanism. Their coherence scales are, however, 
much smaller than those of intergalactic magnetic fields and thus magnetic fields generated from these mechanisms 
could not be directly the origin of large-scale magnetic fields. 

On the other hand, primordial origins, often related with inflation 0, d, [^, H^j [O, [3 o^' second order effects 
through cosmological vector modes [ij, [ij, [l^, [l^ [13, E, 0, have no difficulty in accounting for the length 
of coherence. The observational facts that there exist significant magnetic fields in objects at large redshift may 
support the hypothesis of primordial origin of the large scale magnetic fields [2l|, . If this is the case, it is expected 
that primordial magnetic fields should have formed imprints in the anisotropics of cosmic microwave background 
(CMB) through their stress energy tensor and their Lorentz force on the baryon-photon fluid before cosmological 
recombination. Therefore, it is important to develop the cosmological perturbation theory with primordial magnetic 
fields in order to search for signs of magnetic fields in the observed CMB maps. 

In recent years the effects of stochastic primordial magnetic fields on the evolution of cosmological perturbations 
have been developed independently by several authors (for a revie w I23l). The scalar type perturbations, which is 
related to density fluctuations, have been considered by p^. [25l [26l l27l. l28l [29l. [sot . The vector type perturbations, 
which wouldgive the most dominant contribution to CMB anisotromes at small angular scales, have been studied by 
[31I [3^. [33l. \34. [ssl. [sgI [stI . and also the tensor type perturbations [H, [3^. All of these studies suggest that, from the 
currently available CMB data, the amplitude of primordial stochastic magnetic fields should be at most a few times 
10~^ Gauss or below at the relevant scales (soj . 
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In the present paper we reconsider the scalar type cosmological perturbations with primordial magnetic fields. 
Within the standard cosmological perturbation theory, an initial condition of the perturbation Fourier mode with 
wavenumber k is set when the mode is well outside the horizon {kr ^ 1 with r being conformal time) and when 
the universe is deep in the radiation dominated era, neglecting non-relativistic matter contributions (for example, 
see (ift]). Following this standard practice the initial conditions of density perturbations with primordial magnetic 
fields have been derived [1^, |4il, [42, S, We find, however, that this procedure does not give us a consistent 
initial conditions in the presence of magnetic fields, because of the unwanted cancellation of contributions from the 
magnetic fields and primordial radiations. This cancellation makes the system unstable and violates the constraint of 
perturbed Einstein equations to be satisfied. As we shall show below, this difficulty can be removed by considering 
the significant contributions from the non-relativistic matter at initial conditions. 

The paper is organized as follows. In Section II, we set up basic equations for the perturbation theory inclu ding 
primordial magnetic fields. We adopt the covariant approach to derive equations while equations in Refs. [1^ [s^, Hfl 
are derived in a conventional synchronous gauge. The anisotropic stress and the Lorentz force of the magnetic 
field are defined and the equation of motion for baryons is derived. We also define the spectrum of the magnetic 
field. In Section III, we point out an inconsistency in the previous works and derive the initial condition of the purely 
magnetic mode including the non-relativistic matter contribution. In Section IV we show the numerical calculation 
of CMB and matter power spectra. Finally, we conclude this work. 



II. EQUATIONS 



A. Basic equations 



Here we set up equations. In what follows we take the covariant approach with CDM frame to eliminate the gauge 
freedom ^45, 46, 47]. In this frame, we define variables on the supersurface orthogonal to the CDM 4-velocity u^. Then 
one can define the anisotropic expansion rate (shear) a and the inhomogeneous expansion rate Z from the covariant 
derivative of m^. In the scalar mode, we can neglect vorticity. In addition, we introduce the Weyl tensor, which is the 
traceless part of the Riemann curvature tensor. The Weyl tensor vanishes in the background FRW spacetime. Since 
the magnetic part of the Weyl tensor is negligible in the scalar mode, we define the electric part of the Weyl tensor 
as $. Linearizing the Bianchi identities and Ricci identities, we obtain the following equations for $, a and Z. Three 
propagation equations for <I>, a and Z are: 



1 l~i 1 

^' + n<P + —Kpa^ {ja + RfQf) + —Kpa^{3j - l)i?/7r/ - —npa'Rf^'f = , (1) 

1 



a' + Ha + k^ + -l^Kpa'^RfTTf =0 , (2) 
Z' + HZ + ^Kpa^Rf{Af + 3dPf)^0. (3) 



Two constraint equations are: 



liz-<y) + ^^pa^Rfqf = o, (4) 

1 31~C 
2$ - —Kpa^iRjAf + RfTTf) - —fipa^Rjqj =0. (5) 

The prime " ' " denotes the derivative with respect to the conformal time r and H = a' /a, where a is the scale factor. 
Here we defined 7 as p = (7 — l)p, where p is the total pressure and p is total energy density, and k = 8ttG. The 
subscript " / " means the sum of photon (7), neutrino (i/), baryon (&), CDM (c) and magnetic field (B). The density 
fluctuation Af, heat flux qf and anisotropic stress nf are normalized with their energy density pf except for the 
magnetic field. For the magnetic field variables we normalize As and ttb with the photon energy density pj because 
we consider the case that the magnetic field does not contribute to the background spacetime. Then the total energy 
density p is p = p^ + p^ + pt + Pc- Here we define the energy fraction Rf as R-y = p-y/p, Ru = Pvl P-, RbT = Pb/p, 
RcT = Pc/p and R„i = j{Rb + Rc)- In the deep radiation dominated era, the background energy densities of baryon 
and CDM are neghgible. In this epoch, p ~ p^ + p^, R^ + R^ ~ 1 and all Rfs are constant. Note that Ri, and Rc 
have dimension of (length)"^. Note that in the metric perturbation approach, there are four, not five, equations as 
shown in Ref. [i^. This is because our variables <&, tr and Z are not independent from each other. 

Next we introduce the stochastic magnetic field and evolution equations for each component. In this work, we 
assume that the magnetic field can be treated as the first order quantity and does not contribute to the background 
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evolution. Then the perturbed Einstein equation is described as 

SG^, = 87rG{ST^,+T^^J . (6) 
Here is the energy-momentum tensor for the magnetic field: 



^b'iM = J^, (^S'l - B'{x)Bi(x) ] , (8) 



where Bi{x) is the magnetic field strength at present time and we assume that the conductance of the universe is 
infinite, i.e. Ei = 0. We decompose the space-space part of energy-momentum tensor as 



(9) 



where the first term in the r.h.s. is the trace part and second term is the tracclcss part. The traceless part is the 
anisotropic stress. In the Fourier space, the energy density and anisotropic stress of the magnetic field are defined as 

P^Ab = -Ts°o{k) = <5/Ts^(A:) , (10) 
P^^B = -I (k^y - ^V) Tj,){k) . (11) 

Since T^^^ik) and p-y are proportional to a~'^, and ttb are constant. Then TQ ^j{k) is expressed as 

Tn^jik) = y^jP-yAn - {k% - p^^B • (12) 

If the magnetic field exists, the time evolution of the baryon fiuid is affected by the Lorentz force. The energy 
conservation for baryon is described as 

5T,%,+n\.,,,=o, (13) 

^"^h^n = Pbil'bi + T^Qbi) + ^aUeOTP^i [qu - , (14) 

where Ug is the number density of free electrons and ut is the Thomson cross section. Here we have neglected the 
baryon pressure. Finally, we obtain the equation of motion for baryon in the Fourier space: 



q'b + 'Hqb + aUeaxR iqb- -^Qf j = --kRL , 



L = ^{-Ab + 27Tb) , (16) 

where R = 4p^/(3pf,). The Lorentz force term does not vanish after the recombination because residual free electrons 
and ions still interact with neutral atoms, and thus these particles move together. As wc sec later, the Lorentz force 
have a large effect on the growth of curvature perturbation after the recombination. Under the circumstances where 
the electric field is negligible, the continuity equation for baryon is written in the same manner as the standard one: 

A'b + k{Z + qb) = 0. (17) 

Photons are coupled with baryons through Thomson scattering. The zero and first moments of Boltzmann equation 
for photons are 

A'^ = -k(^Z + q^) , (18) 



q'j = ^(A^ - 277^,) + aneCJT {^qb - q^ ■ 



(19) 
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In the early epoch at which one should impose the initial conditions, the anisotropic stress and more higher multipoles 
of photons are negligible. 

Before the recombination, baryons and photons are tightly coupled, so that qb — 3q-y/4 = v^i,. In the tight-coupling 
epoch, Eqs. (Uni) and are combined to 

<. + ^^.. = |^(A.-3X). (20) 

Other particle species such as neutrino and CDM are treated as coUisionless particles. Since neutrinos are relativistic, 
their evolution should be followed by solving the coUisionless Boltzmann equations: 

^ ~k (^^Z + q,^ , (21) 

ql ^ ^{A, - , (22) 

< = k - ^G(3)) + ^ka , (23) 

G(3)' = hn^ , (24) 

where we set ci,'^ = (Z > 3). CDM can be treated as a nonrelativistic perfect fluid and does not have velocity in 
our frame, qc = 0. Its energy perturbation evolves as 

A'^ + kZ^O . (25) 

B. Spectrum of the magnetic field 

To calculate CMB anisotropics generated from stochastic primordial magnetic fields, we need to specify the spectrum 
of them. It is shown in Refs. [32, 48] that the magnetic field is damped in small scales k > k^, where k^ is the 
wavenumber of damping scale. Here we assume that the magnetic field has power-law spectrum in k < kjj in the 
same manner as in previous works (3^ . Issl . l49l | , 

(S,(k)B*(k')) = {2nf^Ak--S{k-k') , (26) 

with Pij = Sij — kikj, which is the divergence free condition of the magnetic field. There are some ways to define the 
amplitude of the magnetic fields as shown in Ref. [4l[. In this work, we define the amplitude of the magnetic field 
Bx by smoothing at A = IMpc with Gaussian window function in Fourier space, 

Bl = — / d^kAk'"" exp i-X^k^) . (27) 
(27r)-^ J 

Integrating Eq. (I27p . one obtains the two-point correlation of magnetic field, 

(27r)"«+5s2 



i?.(k)i3;(k')) = i2nrP,,-A_L_^_k-^5i^.-k') , (28) 



Since the energy momentum tensor for magnetic field is quadratic in B, one needs to calculate a convolution in order 
to obtain the spectrum of and L. Although many previous works used approximated spectra, Ref. [30j obtained 
the exact spectra for several values of hb- For example, the spectrum for ub = —2.5 is: 



,3! Ai3(fc) |2^_2.5 17fc 



27r2 8007r6 



55 

|i(fc)lL=-2.5=^^|AB(fc)lL = -2.5 ■ (29) 

This spectrum is valid for scales much larger than damping scale, i.e. k ki). Since kjj is sufficiently large, we use 
Eq. throughout this work. 
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-^RjAb 




A, = 


-^R^Ab 





III. DERIVING INITIAL CONDITION 

In previous works [l^, |4l|, \^ initial conditions are derived neglecting the matter contributions. This is a good 
approximation for the adiabatic mode in the standard model. However, in the existence of the magnetic fields, the 
compensation mechanism between the radiation energy perturbation and the magnetic energy density makes the 
matter contributions not to be negligible. We can understand it as follows. The initial conditions for purely magnetic 
mode derived in previous works by neglecting matter contributions are [2^ l4ll 

(30) 
(31) 
(32) 
(33) 

Z = 0{k-'T-'). (34) 

In the radiation dominant epoch, if one substitutes these solutions to Eq. one sees that contributions from 

radiations and magnetic fields compensates each other to obtain 

Z' + nZ+ ^Kpa^RfiAf + 35Pf) = ^{-R„^R^Abt + ©(r^)) + Oik\^) , (35) 

where we neglected the baryon pressure. In the early epoch, r ^ 0, Eq. ((35|) diverges and does not satisfy Eq. ([3]). 
This inconsistency is caused by the matter contribution in r.h.s. of Eq. (j35|) . Thus it is important to derive initial 
conditions including baryon and CDM perturbations. The situation is similar to those of isocurvature models, in 
which the metric perturbations at early times are determined by the non-relativistic matter contributions [50| . 

In order to obtain the appropriate initial conditions with magnetic fields for <I>, a and Z up to the leading order in 
fcr, let us first combine Eqs. ([I])-®. Although the matter contributions make equations complicated, we obtain the 
following three equations in the radiation dominant epoch: a second-order equation for 

(1 + i?™T)$" + (4 + 3i?™T)H$' + l{k^ - + n^)RraT^ * + 0{H)a 

= -^Rf^f [-2H^ + '^n^n' + 2HH" + 2n'^ + {-n^ + 2n'^ + 2H'H")R,nT^ 
+^RfTT'f(2H^ + mn' + {n^ + m.n')RrnT) + -^Rf^'li^ + RmT)H'' 

+ ^(3 + 2R^T)RfTif ~ 2f.R^^TRyAB, (36) 



a third-order equation for a, 



(1 + R,nTy" + (5 + ARrr,T)Ha" + 4H' + 6H^ + (47^2 + AH')RraT + ^ cr' 



3 



+(2n" + mn' + {2h" + m.n')R,nT + y^)^ 

772 

= -rL^{2 + RraT)kRfTTf + —R^rkR^AB 



2 

3 



(^m'^ + mn" + un^n' + [m''^ + mn" + m'^n')Rr,^T^ 2k^f'^f 



- (mn' + m^ + {mn' + m^)R^T) ^Rf^^'f 



-2{l + RraT)n^^RfTr'l, (37) 

and a first-order equation for Z, a and <I>, 

Z' + ^±^HZ + ^-±^k^ ^ ^-±^no . ^^±^i?,., - ^2^^^R,Ab , (38) 

1 ^ RraT 1 + RmT 1 -|- RmT 2k 1 + RmT ^ ^ 2k 1 + RmT ^ ' ^ ^ 
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where we have neglected the baryon pressure and used the adiabatic condition, = Ai, = |Af, = |Ac. As we have 
already pointed out, A^ and irg are constant. The magnetic mode is a particular solution of the linearized Einstein 
equations, while the standard adiabatic and isocurvature modes are the general solutions of them. 

In the purely magnetic mode, photons and neutrinos compensate the energy perturbation of the magnetic field ini- 
tially, i.e. A.^ — = —RjAb + (higher order terms) [41,, i4^, i43i, (4^] , and anisotropic stress of neutrinos compensates 
that of the magnetic field, i.e. tt^ = —^ttb -l-Tr^^^/c^r^, where tt^^^ denotes the coefficient of the ©(fc^r^) term. These 
compensation mechanism and Eqs.([36 |) -([38 |l make it possible to derive $, a and Z up to the second order: 

$ = _ j^R^AskT , (39) 

a = -3i?,7r(2)fcr + ^R^Ask^^ , (40) 



Z = -^R,Ab + (-3i?.7r(2) + i(:^)'i?^AB)fcr , (41) 



where we take into account the matter contribution to the Friedmann equation, namely, a 

K(pbo+Pco) ^2 ^[^Yi being the energy density of species i at present, and 7i ~ r^^ + ^Rm in the radiation dominant 
epoch P,[5l. 

In order to know initial conditions for each component, we need to solve Eqs. p^ - ipS]) . The leading order terms 
of the neutrino perturbation are determined from the compensation mechanism. From Eqs. (j41[) and (j40p . we obtain 
the neutrino perturbation up to the second order: 



A^ = -R^Ab + -RrnR-yABT , (42) 



= (^R^Ab - S^T^B^ kr + ^R„,R^ABkT^ , (43) 

= _ ( IrA^ _ -^TTB + -R.TT^'A -k'r^ , (44) 

i?^ V15 21 i?,. 5 y 2 ' ^ ' 

Gi'^ = ~\^^-Bkr . (45) 

Again, the definition of tt^^) is tt^, = -^t^b + Tr^^^fc^r^. Then Eq. leads to 

(2) 1 R^ IARi^Ab - 557r_B , . 

~ 42i?^ 4i?^ + 5 ■ ^ ' 

The equation of motion for photon-baryon fluid, Eq. (I20p . imply that 

\{RuAb - 2nB)kT + q^^'hr^ . (47) 
From Eqs. ([4]) and (I4ip and tight coupling approximation, qb ~ Sq-y/A, we can obtain q^^: 

-,(2) 

= ^R„,R^Ab - ^^{R,Ab - 27rs). (48) 
The Boltzmann equations and the solutions of Z, q^ and q^ give A^, A^ and Ac up to the third order: 

A, = -R^Ab + ^i?™i?^ABr + (-2^(2)^^ _ + \{R^Ab - 2^i:B))k\^ , (49) 



q, 3, 



Rfqf = [ -RrnRuR-yAB + R^q)y^ + -Rb{R^AB - 2tTb) ) kr^ = -R,nR-yABkT' 



^Ab = -R-,Ab + ^RmRyABT + (^-27:^^^ R, - ^(^^y R^Ab - ^{R^Ab - 2nB))pT^ , (50) 



Ac = ^Ir,Ab + ii?„,i?,Asr + (-^vr(2)i?, - ^ (:^) 'i^^A^) fcV^ 



(51) 
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In the numerical calculation, we use the curvature perturbation r], i] = ~(2^ + (t' /k), in place of $. Its initial condition 
is 

r; = -6i?.^(2) + i:^i?,ABfcr-l(:^)\AsfcV . (52) 

Since this magnetic mode is a particular solution of the linearized Einstein equations, the general solution can be 
the sum of standard adiabatic mode and purely magnetic mode. In such cases, the total temperature perturbation 
A'°* is represented as 

A*°' = A^'^' + A^ . (53) 

Here A'^'*' is a temperature perturbation from the adiabatic mode, which is calculated with standard adiabatic initial 
conditions (40j . and A^ is from the purely magnetic mode. Then the ensemble average is 

(A'°*A*°**) (A'^'^'A^'^'*) + (A^A^*) + (A'^^'A^*) + (A^A^"^'*) , (54) 

where the latter two terms are the correlation between the adiabatic and the purely magnetic modes. In what follows 
we study the three cases, namely, fully correlated case, anti correlated case and uncorrelated case. 



IV. RESULT AND DISCUSSION 



A. Matter Contribution 



The new initial conditions derived in the above are totally different from those used in the previous works. Since 
we do not neglect matter contribution, A's and rj have terms of 0{RmT). To see the importance of RmT terms, we 
calculate the total energy perturbation up to the second order: 



RfAf = --R^R^Abt + ( ~2R,^^'^k' + —R'^R^Ab ) t' (55) 
^R^Ab + ^R^R^Abt'^R^t - 2i?,^(2)^2^2 _ (5g) 

We define a wavenumber /cmat as the scale in which absolute values of two factors in 0{t^) term in Eq. ((55|) are equal 
to each other, 



2 

mat 



—RI^RjAb 



(57) 



then we obtain /cmat — l-2i?m. If A; < fcmat, Rf^f is dominated by the matter contribution. Since the wave number 
^mat correspond to the multipole Zmat ~ ^mat'na ~ 80, where tq is the present conformal time, we cannot neglect the 
matter contribution in lower multipoles. When k > k^-nnti we define as the conformal time when the absolute 
values of two terms in r.h.s. of Eq.(|56p are equal to each other, i.e. 

kr ^RmR^Ask f2— + 175:^r' f58) 

" - 36 I 77(2) I ^^fc2 + ml^R^AB -\ R^^ k ) ' ^^^^ 

where we assumed 1 3> A:t > k^^^^r ^ RmT- When kr is smaller than kr^, RmT term in Eq. (|56p . which comes from 
the matter contributions, can not be negligible and should be incorporated appropriately. In the radiation dominant 
epoch, R„i is represented as 

3 Qbh^+n^h^ 100 km sec-i Mpc-i 

" -^V^^h^ + n^h^ ^""P^" ^ 

^ 0.02 + 0.11 V 4.3 X 10-5 / ^ ' 

Then we obtain kr^ as shown in Fig. [1] In small scales, fcr,„ become smaller and we can neglect the matter 
contribution. However, if we start to integrate equations from the conformal time much smaller than r™, the matter 
contribution plays an important role in the magnetic mode. For example, if one sets initial conditions at fcr ~ 0.01, 
then from Fig. [T]we find that the modes k < 0.3 suffer from the matter contributions. These modes correspond to 
the angular scale £ < 4000, at which the difference should be significant as shown in Fig. [T] 



8 



B. Numerical Calculation 

The equations and initial conditions derived in the previous sections can be used to calculate CMB and mat- 
ter power spectra numerically. We calculated them with accordingly modified CAMB code [5l|. In all of 
our calculation we fixed cosmological parameters to the best-fitting values to the WMAP-5yr data [s^l, namely 
{u;b,uJc,h,Tc,A%ns) = (0.0227,0.1099,0.719,0.087,2.41 x 10-^,0.963), where ujb and uJc are the energy densities of 
baryon and CDM, respectively, h is the Hubble parameter, Tc is the optical depth, A|j and ris are the amplitude and 
the spectral index of primordial curvature fluctuations, respectively. 

In Fig. [U we compare the CMB spectrum of purely magnetic mode calculated with our new initial conditions 
and that with previous initial conditions, in which the matter contributions were omitted. At higher multipoles, 
two spectra converge with each other asymptotically. This is because the matter contributions are negligible in very 
small scales as shown in the previous subsection. However, at large angular scales, we cannot neglect the matter 
contribution and the correct initial condition yields larger amplitude. 

The perturbed Einstein equation gives four equations and two out of four are independent. In order to check 
the consistency of our numerical calculation, we picked up six kinds of different pairs of two equations from the 
four independent equations and observed that all of the results coincide with each other. Note that the old initial 
conditions do not satisfy the Einstein equation. We found that different pairs of perturbed Einstein equations yield 
different results if we start with the old initial conditions. 

In Fig. [21 we plotted the baryon heat flux qb and curvature perturbation 77 normalized by the square root of the 
power spectrum, -^Z | Ab p, qb = qb/y/ \ Ab p and i) = rj/y/ \ Ab p. The characters of the growth of perturbations are 
different between k > fcrec and k < k-ceci where k-^cc is a wavenumber which crosses the horizon at the recombination 
epoch. In the case of fc > fcrec, the perturbation enters the horizon before the recombination. After the horizon- 
crossing, the growth of the baryon velocity is suppressed and shows oscillatory behavior by the photons pressure 
through Thomson scatterings. However, after the recombination, qb grows suddenly by the Lorentz force. This qb 
evolution enhances 77 because the source of the curvature perturbation is the total velocity field, i.e. 

77' = -^Kpa^Rfqj . (60) 

On the other hand, at scales where the waves enter the horizon after recombination, k < k^cc, the baryon velocity does 
not undergo the suppression and there is no sudden growth of potentials (dash-dotted (blue) line in the left panel of 
Fig. [2|). 

In Fig. [3l two metric perturbations in conformal Newtonian gauge are plotted. Again, we plotted normalized 
potential, and "0, with respect to the square root of the power spectrum. Potentials and "0 are defined as 

ds^ = a^(-(l + 2^j)dT^ + (1 - 2(j3)dx^) , (61) 

which is the same definition as in Ref. [40]. The relations with the variables used in this paper and their initial 
conditions are given as 

= 77 - Ina = -3i?,7r(2) + ^R^AekT , 

^ ^ hfj' ^jia) = ^Qjl^T^i^) + ^^jl ABkT . (62) 
k 8k 

In the absence of the magnetic field, two potentials decay to zero, once the perturbation enter the horizon in radiation 
dominant epoch (thick lines in Figs. [5] and [3]), and are constant in time in matter dominated epoch. However, in the 
purely magnetic mode, potentials grow after the recombination. 

Next we study the growth of curvature perturbations in the adiabatic mode correlated with the magnetic mode. In 
Fig. m the growths of 77 in full and anti correlated cases are plotted. If the adiabatic mode is fully correlated with 
the magnetic mode, 77 grows after the recombination for k > fcrcc, while 77 decays in the anti correlated case. This 
growth of the curvature perturbation is directly seen in the CMB and matter power spectra as shown in Fig. [S] For 
the CMB spectrum (three panels in Fig. [5]) , primary standard adiabatic mode and uncorrelated mode have a similar 
feature in shape and almost degenerated. Therefore, if the adiabatic mode is fully correlated with the magnetic mode, 
the gravitational potential becomes deeper and the amplitudes of spectra are increased. On the other hand, the anti 
correlated magnetic mode decays the potential and amplitudes of spectra become smaller. 

For the matter power spectrum, on the other hand, the difference shows up at small scales because the Lorentz 
force from magnetic fields newly induces the density perturbations dominantly at small scales after cosmological 
recombination. For both correlation cases with 300 nC magnetic field, the linear power at fc > IMpc^^ is dominated 
by the perturbations induced by magnetic fields. 
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FIG. 1: Left: The conformal time at which the matter contribution becomes significant, fcT,„. If kr < krm, the matter 
contribution can not be neglected. Right: Comparison of the CMB power spectra derived from the initial condition with 
matter contribution and from the one without. The new initial condition leads the larger amplitude. Magnetic field parameters 
are Bx = InG and ub = —2.5. 
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FIG. 2: The normalized baryon heat fiux qt (left) and curvature perturbation ?) (right) of the purely magnetic mode. The 
baryon velocity evolves by the Lorentz force after recombination. This enhances the curvature potential at small scales. Each 
thin line represents the purely magnetic mode with different scale: solid (black), dashed (red), dotted (green), dash-dotted 
(blue), double dotted (magenta), and dashed-double dotted (light blue) lines corresponds to the modes with k = 10, 1, 10~^, 
10~^, 10""^, and 10~* Mpc~^, respectively. In the right panel, a thick solid line represents the standard adiabatic mode with 
k = lOMpc"^, which decays to zero in the radiation dominated era and stays constant in the matter dominated era. 



V. CONCLUSION 



In this work, we derived the initial conditions for density perturbations in the existence of the primordial magnetic 
field. Because the compensation mechanism between radiation perturbation and the magnetic field makes it impossible 
to neglect matter contributions even in the early universe, we derived initial conditions including matter contributions 
to the density field. The initial condition derived in this paper fully satisfies the linearized Einstein equations, and 
thus enables us to solve the system numerically in a stable and consistent manner. Then CMB and matter power 
spectra are calculated, and the evolutions of perturbations are presented in detail. We found that it gives the larger 
amplitude of the CMB angular power spectrum by at most a order of magnitude at large scales I < 4000, compared 
with the initial condition in the literature. 

In the purely magnetic mode, the potentials grow suddenly after recombination for k > k^cc- This effect enhance 
the amplitude of the matter power spectrum if the adiabatic and magnetic mode are fully correlated, and decreases 
the amplitude in the anti correlated case at intermediate scales. For much smaller scales A; > 1 Mpc~^, the spectrum 
is dominated by the magnetic mode. 
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FIG. 3: Gravitational potentials <j) and ip of the purely magnetic mode in conformal Newtonian gauge. Potentials grows 
suddenly at small scales after the recombination. Each thin line represents the purely magnetic mode with different scale: solid 
(black), dashed (red), dotted (green), dash-dotted (blue), double dotted (magenta), and dash-double dotted (light blue) Unes 
corresponds to the mode with k = 10, 1, 10~^, 10~^, 10""^, and 10"* Mpc~^, respectively. For both panels, thick solid lines 
are the standard adiabatic modes with k = lOMpc"'^. 




FIG. 4: The growths of r) in fully (left) and anti (right) correlated cases. The magnetic field enhances the curvature perturbation 
in the fully correlated case. On the other hand, anti correlated magnetic field diminishes ry. Each thin line represents different 
scale: solid (black), dashed (red), dotted (green), dash-dotted (blue), double dotted (magenta), and dash-double dotted (light 
blue) lines corresponds to the mode with k — 10, 1, lO"'^, 10~^, 10~^, and lO"** Mpc~^, respectively. Parameters for magnetic 
fields are Bx = 300nG and ub ~ —2.5. 
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